Epidemic spreading in evolving networks 
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A model for epidemic spreading on rewiring networks is introduced and analyzed for the case of 
scale free steady state networks. It is found that contrary to what one would have naively expected, 
the rewiring process typically tends to suppress epidemic spreading. In particular it is found that as 
in static networks, rewiring networks with degree distribution exponent 7 > 3 exhibit a threshold 
in the infection rate below which epidemics die out in the steady state. However the threshold is 
higher in the rewiring case. For 2 < 7 < 3 no such threshold exists, but for small infection rate the 
steady state density of infected nodes (prevalence) is smaller for rewiring networks. 
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I. INTRODUCTION 

Epidemic spreading can be thought of as occurring on 
complex networks where the nodes of the network repre- 
sent individuals and the links represent various interac- 
tions among those individuals. For example the spread- 
ing of diseases can be thought of as occurring over the 
network of human contacts |l| and the spreading of com- 
puter viruses as occurring over the internet Models 
of epidemic spreading over networks have been studied 
extensively in recent years (for reviews see [3, [B[ ) . Typ- 
ically, the underlying network in these models is consid- 
ered to be static while the state of the individuals resid- 
ing on its nodes can change from infected to non-infected 
according to some dynamical rules. One is then inter- 
ested in studying the evolution of an infected region in 
time, the average density of infected nodes in steady state 
(prevalence) and the way they are affected by the statis- 
tical properties of the network and the infection rates. 

In general, networks can be characterized by the con- 
nectivity of their nodes. The connectivity (degree) k of 
a node is defined as the number of links connected to 
the node. The degree distribution of a network V{k) 
is defined as the probability of a randomly chosen node 
to have a degree k. Many networks such as social net- 
works, the internet and the World Wide Web (WWW) 
have been found to be scale free (SF) meaning that 

the degree distribution follows a power law 



V{k) ~ . 



(1) 



In the thermodynamic limit one can divide SF networks 
into two classes based on the exponent 7. For 7 > 3 the 
second moment of the degree distribution is finite and as 
such the system exhibits finite degree fluctuations. For 
2 < 7 < 3 the second moment diverges resulting in in- 
finitely large degree fluctuations. In the present study we 
only consider networks with a finite degree distribution 
corresponding to 7 > 2. Interestingly, many real net- 



works have been measured to belong to the second class 
having 2 < 7 < 3 0. 

Studies of models of epidemic spreading over static net- 
works have shown that in networks for which 7 > 3, 
the prevalence, p vanishes for sufficiently small infection 
rates A. The prevalence become non-zero only beyond a 
threshold rate A c . On the other hand for networks with 
2 < 7 < 3, for which the second moment of the degree 
distribution diverges, the prevalence is non-zero for any 
infection rate, and no threshold exists fl0l - [l3l ]. Thus, epi- 
demics are easier to stop in static networks with 7 > 3. 

In many cases networks are not static but rather evolve 
in time, for example via rewiring processes. Steady states 
of rewiring networks have been studied in the past. It has 
been shown that depending on the average degree and the 
rewiring rates, networks may reach an SF steady state, 
with an exponent 7 which can be expressed in terms of 
the dynamical rates [T^ - [l6| . 

In the present paper we consider epidemic spreading 
over rewiring networks. On such networks, the disease 
can spread at a given time through the links which are 
present at that time. We find that as in the static case 
a non- vanishing threshold value of the infection rate, A c , 
exists for 7 > 3. Below this threshold the prevalence 
(fraction of infected individuals) vanishes while above it 
the prevalence is non-zero. For 2 < 7 < 3 no such thresh- 
old exists and the steady state prevalence is non-zero for 
any A > 0. However, contrary to what one would have 
naively expected, epidemic spreading in our model is not 
necessarily enhanced by the dynamics of the network. 
For 7 > 3 the threshold A c is found to be larger than that 
of the corresponding static network. Also, for 2 < 7 < 3 
the prevalence at small A is found to be smaller than that 
of the corresponding static network. 

The paper is organized as follows: in section [TT] we re- 
view known results on epidemic spreading in static net- 
works and on networks with rewiring dynamics. In sec- 
tion IIIII we study epidemic spreading on evolving net- 
works using mean field calculations and numerical simu- 
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lations. Our results arc summarized in section HV 



II. REVIEW OF KNOWN RESULTS 
A. Epidemic spreading in static networks 

A number of models of disease spreading have been in- 
troduced and studied in the past. In the present work 
we use the Susceptible Infected Susceptible (SIS) model 
0, !, H-GH, [3 [l3- In this model a healthy individ- 
ual, with respect to the disease, may be infected through 
interaction with diseased individuals. Meaning, that a 
susceptible node may be infected through a link connect- 
ing it to an infected node, which we will refer to as his 
neighbor. Once an individual is infected he may become 
susceptible again by being spontaneously cured from the 
disease. The curing process does not immune the indi- 
vidual and it can be reinfected. 

The continuous time dynamics of an epidemic in the 
SIS model is defined by two stochastic processes using 
two parameters: 

A- Infection rate 
S - Rate of recovery 

An infected node is spontaneously cured with a rate S 
which we choose to be equal to 1 by adjusting the time 
scale. On the other hand a susceptible node gets infected 
with rate A from each of its infected neighbors. Thus, the 
rate a node is infected depends linearly on the number 
of infected neighbors. This model of infection is different 
from the model explored in EH EH, E3 where the 
infection rate is independent of the number of infected 
neighbors. However, both models behave similarly near 
the threshold for an endemic state and we expect our 
conclusions to hold for both models. 

The problem is addressed using a mean field (MF) ap- 
proach and numerical simulations. The MF approach ne- 
glects correlation in infection between nodes in the sense 
that for any pair of nodes i,j we have (f]iT)j) = {rji){rij) 
where r] = 0, 1 is a parameter indicating whether a node is 
susceptible or infected, respectively. As an order param- 
eter we use the prevalence of the disease, the density of 
infected nodes in the network, defined as p = Ni n { ecte ^/N. 
Hence, our problem is reduced to a contact equation for 
the order parameter p. Since we are interested in for- 
mulating the problem for any degree distribution, as was 
previously done in [ToL ET| , we shall distinguish between 
nodes of different degree by defining p/~ as the fraction of 
diseased nodes of degree k. The total prevalence is thus 
given by 




u(k)v(k' 



FIG. 1: Rewiring of a link from a node with degree k to a 
node with degree k' with a rate u(k)v(k'). 



With t being the density of "infected links" defined as 



Note that 



E 

fc=0 



V*(k) 



kV{k) 



-Pk 



kV(k) 



(4) 



(5) 



is the degree distribution of a randomly chosen neighbor 
of nodes. Thus (j4j) gives the probability that a randomly 
chosen end of a randomly chosen link is infected. In the 
steady state, a non- vanishing solution for the prevalence 
is possible only for infection rates greater than (see [l(| 

m 



A, 



{kl 



(6) 



For infection rates above the threshold, A > A c , a finite 
fraction of the nodes is infected while for A < A c the dis- 
ease dies out 1 . 

For Erdos-Renyi (ER) networks, which obey a Poisson 
degree distribution [l8j, the threshold can be rewritten 
in the form A c = l/((fc) + 1). Moreover, for SF networks 
with 2 < 7 < 3 the second moment diverges and as a 
result the threshold vanishes. As a consequence such a 
system will always reach an endemic steady state for any 
non zero infection rate A > 0. 



B. networks under rewiring dynamics 

During rewiring dynamics of a network the number 
of nodes and the number of links are unchanged but the 
links are stochastically detached from one node and reat- 
tached to another. In our model the process of rewiring 



fc=0 



The MF contact equation has the following form 

dpk 
dt 



-pk + \kl(l - p k ) 



(2) 



(3) 



1 In the thermodynamic limit this is a transition to an absorbing 
state, but for finite size systems the only true steady state is one 



with zero prevalence 1121 . As a result, in finite networks there 
is no true threshold but a crossover infection rate which can be 
calculated for a quasi-stationary state. 
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a randomly chosen end of a link from a node with degree 
Ho a node with degree kl occurs with rate u{k)v{k'). A 
schematic representation of the process is given in Fig. [TJ 
These rates determine the steady state degree distribu- 
tion of the network through the relation 



n^(fc') 



V(k) 



n£=i«(*o 



V(0) 



(7) 



which can be derived from the master equation for the 
node degree distribution [Tif . 

Under such rewiring dynamics, the resulting networks 
are uncorrelated in the sense that the joint probability 
n(k, k') that the ends of a randomly chosen link are nodes 
of degree k and k 1 factorizes to 



7r(k,k') = V*(k)V*{k') 



(8) 



By choosing the proper attachment and detachment 
rates one can create an evolving network with a constant 
size and any desired degree distribution. One such choice 
of rewiring yields an evolving ER type network. This is 
achieved by choosing a link at random and rewiring one 
of it's randomly chosen ends to a randomly chosen node. 
This rewiring scheme has a constant attachment rate and 
a linearly preferential detachment rate. One can easily 
verify by using ([7]) that the choice 



1 

N 

u(k) = k 



v(k) 



(9) 



indeed yields a Poisson degree distribution. 

Through the use of such rewiring dynamics we can cre- 
ate an uncorrelated SF network with any desired expo- 
nent in the power law distribution. In what follows we 
work with rewiring dynamics similar to that of zero-range 
processes (ZRP) [la, Gil > wnere the rewiring rate does not 
depend on the destination site, i.e., v(k) — 1/N. As a 
further specification we consider detachment rates of the 
form 



u{k) = l + r 
k 



(10) 



with b as a parameter of the dynamics. In this case, for 
a specific choice of the average number of links (given 
by (k) = 1/(6 — 2)) the underlying zero-range process 
exhibits critical behavior in which the steady state degree 
distribution is a power law V(k) ~ k~ b at large k. At 
lower average link numbers the steady state distribution 
decays exponentially with k while at larger averages a 
hub becomes present which is linked to a finite fraction 
of the nodes in the network [l5|, [l|| . 

In order to be able to control the critical value of (k) 
for a given value of b one can make a slight modification 
in the dynamics by considering 



u(jfe) 




< k < k Q 
k > k 



(11) 



Due to the same asymptotic behavior of (|T0|) and (fTTj) 
this modification does not change the power law tail of 
the stationary degree distribution. In this case the crit- 
ical value of (k) can be obtained numerically Note 
that for rewiring rates of the form (fTOj) and (fTTj) , (u) = 1 
at criticality [l5| . 

It is important to note that the dynamics, as defined, 
allows for multiple link between two nodes (melons) and 
links that connect a node with itself (tadpoles). By not 
allowing for melons and tadpoles we are introducing an 
effective preferential attachment rate v(k) = 1 — k/N 
as opposed to a constant rate as given in (fTU|) . This 
rate takes into account the fact that the neighbors of 
a node of degree k are not available as target nodes for 
the rewired link. The preferential attachment rate means 
that a highly connected node has a lower rate of attach- 
ment than a node with a lower connectivity and induces 
disassortative, or negative correlations. It can be shown 
using (|7|) that this attachment rate imposes a Gaussian 
cutoff on the degree distribution of the form 



V\k) 



T V(k) 



(12) 



where V(k) is the degree distribution for similar dynam- 
ics which allows for tadpoles and melons. For ER type 
networks and for SF networks with 7 > 3 the fraction of 
melons and tadpoles vanishes in the thermodynamic limit 
[l9l | . However, for SF networks with 7 < 3 the number of 
melons and tadpoles diverges and cannot be neglected. 
Since the infection process was taken to depend linearly 
on the number of infected links, the problem could be re- 
stated by choosing networks with weighted links and an 
infection process which depends linearly on the weight of 
the link. 



III. EPIDEMIC SPREADING ON EVOLVING 
NETWORKS 



Our aim is to consider a model of epidemic spreading 
on a network which is changing in time. As a conse- 
quence, a given node is no longer connected to a static set 
of neighbors but to a dynamic one, and the degree k of the 
node also fluctuates. Previous work on epidemic spread- 
ing on evolving networks poU^ concentrated mainly on 
rewiring dynamics resulting from the adaptation of the 
network to the disease. In these models the rewiring dy- 
namics depend on the state of the nodes, i.e. the infection 
process. We consider models where the rewiring dynam- 
ics is independent of the infection process. To be more 
specific, we consider a ZRP-likc rewiring dynamics for 
the network with rates of the form pT|) and v(k) = 1/N. 
For a specific value of (k) this results in a power-law 
degree distribution for the steady state of the network 
(which depends on fc ) with V(k) ~ k~ b . In addition, 
we introduce a parameter v, which describes the overall 
timescalc of the rewiring process as compared to that of 
the infection process. The rewiring rate from a node of 
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degree fc then becomes vu(k). For v — the model re- 
duces to epidemic spreading on a static network, whereas 
for v — > oo, due to the fast mixing we expect a mean-field 
like behavior for the infection process, where neighbors 
change very rapidly. We note here that we always con- 
sider the case where the network is in a stationary state 
with respect to the rewiring dynamics, which requires a 
diverging equilibration time for v — > 0. 



A. Mean-field results 

In order to account for the rewiring dynamics ([3]) has 
to be modified as follows: 

= -Pk + \k£(l - p k ) - vpk{uk + (u)) 



+ v 



V(k + 1) V{k-l). , 



(13) 



By multiplying (|13[) with V(k) and summing up for all k 
one obtains 



p = \(k)t{l-l) 



(14) 



in the stationary state. For infinitesimal p and I (at the 
threshold) this reduces to 

p = K{k)L (15) 

Note that one can rewrite definition @ of i as 

e = p(k) in{ /(k), (16) 

where (-)mf denotes an average in the ensemble of in- 
fected nodes. We define the average of a quantity x in 
the ensemble of infected nodes as 



(x} in i = P 1 ^2xp k V(k). 



Using this, (fl~5j) takes the following simple form 



(k) 



inf- 



(17) 



(18) 



On the other hand, multiplying (|13[) by kV{k) and 
summing up over all k one obtains using (|15|) for the 
steady state 



A,- 1 = + ^((u}-( U ) i „ f ). 



(19) 



Whether or not there exists a non-zero infection rate 
threshold can be easily deduced from this equation. For 
the rewiring rates (jlip both (u) and (u)- m f are finite. 
Thus there exists a finite positive threshold as long as 
(k 2 ) is finite, namely for networks with 7 > 3. In this case 
the rewiring rate v affects the threshold quite strongly. 
On the other hand for networks with 2 < 7 < 3, (fc 2 ) 
diverges and A c vanishes. 



It is obvious that for v — > (jl9|) reduces to A c = 
(k) I '(fc 2 } as discussed before. In the other extreme case, 
when v — ¥ 00, due to the fast rewiring we expect that 
the degree distribution of infected and non-infected nodes 
become identical. This would imply that p k = p for 
all fc and (fc)i n f = (fc). Therefore, based on (|T5)) . one 
has A c = (k) . Note also that in this infinite rewiring 
limit the rhs of (|19|) is nontrivial, since v — > 00, while 
(u) — (u}i n f is expected to vanish. 

It is important to note that whereas for small values 
of v the MF approximation, that we use throughout this 
section, is not necessarily valid. However in the v — 5- 
00 limit the fast rewiring ruins all the correlations in 
the system, and the MF approximation is expected to 
become asymptotically exact. 

As shown by (|18[) the threshold is determined by the 
degree distribution of infected nodes at the transition 
point. In the following we attempt to get a deeper un- 
derstanding of how this distribution changes with the 
rewiring rate v. For this reason we define 



,. Pk 

r k = hm — , 



(20) 



and assume that this quantity is finite for all fc. This 
implies the following normalization for r^: 



J2V(k)r k = l. 



k=0 



It is easy to see that 



V(k)r k = V{k) 



(21) 



(22) 



is the degree distribution of infected nodes close to the 
threshold. 

Equation Q13[) together with the steady state relation 



V(k) _ Mfc+i 
V(k + 1) ~ (u) 



(23) 



implies 



= -p k + Xk£(l - p k ) 

+ v [(u) Pk+1 + u(k)p k -i - ((«> + u{k))p k ] (24) 

for the steady state. Inserting ([20]) into (|24| and using 
(f2Tt]) and (|T5|l one obtains the following set of equations 
for r k at the transition point: 



= jj-j - r k + v [r k+ i(u) + r k -iu k - r k (u k + (u))} . 

(25) 

After solving the above set of equations for r k one can 
determine A c from (fT8|) and (f20l as 



X- 1 = J2kr k V{k) 



(26) 



One can immediately see that in the v — ¥ 00 limit the 
solution of (|2"5"|) is r k = 1. This implies V = V which 
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FIG. 2: (Top) rfc is plotted against k for various rewiring 
rates. The set of equations (|25|) is solved numerically with 
b = 3.5 and ka = 15 for the critical case, which corresponds to 
(u) = 1, (k) = 4.917454, and (A: 2 ) = 99.39410. (Bottom) The 
threshold A c as a function of v. It was calculated numerically 
using equation (|26|) with the same parameters. 



results in the already noted limiting behavior with A c = 
l/(fc). On the other hand, in the case of a static network, 
where v = 0, pk is proportional to k near the threshold, 
implying = k/(k) and V = V*, which results in A c = 
(k)/{k 2 )- 

Numerical solutions of (|25p for intermediate values of v 
are shown in Fig.[5J One can see that for a large but finite 
v there are two crossover values of k. For k below some 
value k\ one has Tk = 1 and V = V (infinite rewiring), 
whereas for k larger than some other value &2, one has 
fk = k/(k) and V = V* (static). Between k\ and ki 
we find an intermediate regime, which connects the two 
extreme cases. The crossover values k\ and hi increase 
with increasing v. 

Using these numerical solutions for r& we calculated 
numerically the effect of rewiring on the threshold A c by 
using (|26|) . Results are shown in Fig. [2] One can clearly 
see that as the rewiring rate increases the threshold in- 
creases. 



B. Simulations 

As discussed above, for SF networks with 2 < 7 < 3 
there is no threshold in the infection rate, and the preva- 
lence is non-zero for any A > 0. For SF networks with 
7 > 3 a threshold exists such that for an infection rate 
below A < A c the prevalence is zero. The prevalence 
corresponding to such networks is studied in this section 
using numerical simulations of finite networks. 

The networks were constructed using two sets of pa- 
rameter values for the rewiring dynamics As an 
example of a network with 2 < 7 < 3 we use the dy- 
namics with v(k) = 1/N, u(k) = 1 + b/k (ko = 1) and 
b = 2.5 . As described previously, (k) = 1/(6 — 2) = 2 
corresponds to the critical average value of the underly- 
ing zero-range process for which the steady state degree 
distribution is a power law with 7 = 2.5. The resulting 
prevalence as a function of the infection rate is plotted 
in Fig. G3a) for various rewiring rates. 

As an example of a network with 7 > 3 we used the 
dynamics JTTJ with (k) « 4.917, b = 3.5 and k a = 15, 
corresponding to the critical average value of the under- 
lying zero-range process for which the steady state degree 
distribution is a power law with 7 = 3.5. The resulting 
prevalence as a function of the infection rate is plotted 
in Fig. EJa) for various rewiring rates. 

For a network of finite size, there is no true thresh- 
old but a crossover infection rate which is obtained for 
a quasi-stationary state. In simulations one can identify 
this crossover value X C (N, L) as the point where the (nu- 
merically obtained) derivative dp/dX takes its maximum. 
With such a definition \ C (N,L) — > A c in the thermody- 
namic limit. One can see in Fig. |3|b), corresponding to 
2 < 7 < 3, that as the rewiring rate increases X C (N,L) 
increases from approximately (k) / (k 2 ) w 0.04 towards 
l/(k) = 0.5. Similarly, one can see in Fig. HJb), corre- 
sponding to 7 > 3, that as the rewiring rate increases 
X C (N,L) increases from approximately (k) / (k 2 ) 0.07 
towards l/(fc) = 0.2. 

In the simulations a very weak external source of infec- 
tion was introduced in order to prevent the system from 
fluctuating into the absorbing state. There are several 
other methods of simulating an absorbing phase transi- 
tion and computing from it the value of the threshold 
which are reviewed in {251 ]. 

Note that for the considered value of b = 2.5, (k 2 ) di- 
verges in the thermodynamic limit, therefore, based on 
the MF results, the threshold is expected to vanish for 
v = 0. However, for a finite system one expects a finite 
crossover value for A, which is of order l/(k 2 ). On the 
other hand, for v (k 2 ) the crossover value should in- 
crease up to l/(fc). It is interesting to examine how the 
crossover A changes if the rewiring rate is of order (k 2 ). 
To this end we performed simulations with v ~ (k 2 ). We 
found that the threshold scales as Aj 1 ~ (k 2 ). The cor- 
responding data collapse is presented in Fig. [5] where the 
prevalence is plotted as a function of a scaled infection 
rate X c (k 2 } for networks of different sizes with a rewiring 
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FIG. 3: The prevalence p (a) and the derivative of the preva- 
lence with respect to the infection rate dp/dX (b) are plot- 
ted as a function of the infection rate for a network of size 
N — 200, (fc) — 2 and rewiring dynamics with b — 2.5 
and fco = 1 for which (fc 2 ) « 58. Curves for rewiring rates 
u = 0, 10, 10 2 , 10 3 , 10 4 , 10 s are represented by full lines from 
left to right respectively. For clarity, the derivative was cal- 
culated using smoothed interpolated data of the prevalence. 
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FIG. 5: The prevalence as a function of the scaled infection 
rate X c (k 2 } for a network of different sizes N, (k) = 2 and 
rewiring dynamics with b — 2.5 and ko — 1 plotted for 
different rewiring rates where for each network of size N the 
rewiring rate is v = {k 2 }. The behavior near the threshold is 
replotted on a finer scale and is given in the inset. 






FIG. 4: The prevalence p (a) and the derivative of the preva- 
lence with respect to the infection rate dp/dX (b) are plot- 
ted as a function of the infection rate for a network of size 
N = 200, (k) = 4.9 and rewiring dynamics (TTTj) with b = 3.5 
and fco = 15 for which (fc 2 ) w 73. Curves for rewiring rates 
v = 0,10, 10 2 ,10 3 are represented by full lines from left to 
right respectively. For clarity, the derivative was calculated 
using smoothed interpolated data of the prevalence. 
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FIG. 6: The prevalence as a function of the infection rate for 
a network of size N = 200, (fc) = 2 and rewiring dynamics 
pip with b = 2.5 and fco = 1. The simulation results for 
v = and v = 10* are compared with the numerical solution 
of (0. 



rate equal to the second moment of the degree distribu- 
tion v = (k 2 ). 

In Fig. |5]and Fig.[7]the simulation results are compared 
to the numerical solution of (fl"4| for b = 2.5 and b = 3.5 
for both a static network [y = 0) and for v = 10 4 . The 
numerical calculation of the MF contact equation was 
carried out by solving (|13[) for each infection rate. The 
degree distribution used in the calculation was taken from 
the simulation results. For both the b = 2.5 and b = 
3.5 cases the MF solution for v = 10 4 agrees quite well 
with the simulation results which supports our general 
argument that as the rewiring rate increases, compared 
to the infection process, both the degree of a node and 
fluctuations average out such that the MF approximation 
better describes the process. 



= MF v-0 
q g _ • simulation v=0 

MFv=10 4 
nc. ♦ simulation v=10 4 



FIG. 7: The prevalence as a function of the infection rate for a 
network of size N = 200, (fc) « 4.917 and rewiring dynamics 
(|11[) with b — 3.5 and fco = 15. The simulation results for 
v = and v = 10 4 are compared with the numerical solution 

of (rjU). 
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IV. CONCLUSIONS 

The effect of network dynamics on epidemic spreading 
has been studied using mean field analysis and numer- 
ical simulations. In particular we considered epidemic 
spreading over SF networks with rewiring dynamics. 

We have shown that the introduction of rewiring affects 
the threshold for an endemic state of a network. This is 
a surprising result that an evolving network is fitter with 
respect to disease the faster it is rewired. This result is 
general to any network with a general degree distribution. 

One can understand this counter intuitive result by as- 
sociating the second moment of the degree distribution 
with the heterogeneity of a network. The more heteroge- 
neous is a network the larger is the fraction of highly con- 
nected nodes which mediate the infection process. The 
introduction of rewiring effectively averages out the het- 
erogeneity and creates an effective homogeneous network, 
with respect to the infection process, where each node has 
an effective average degree k = (k). 

Different networks differ in the rate of rewiring that is 
required for a change of the threshold. We have shown 
that for networks with different degree distributions the 



relevant quantity is the second moment of the degree dis- 
tribution. Only if the rewiring is larger than the second 
moment A c > (k 2 ) then the threshold is affected and is 
increased from A,: 1 = (k) / (k 2 ) to A^ 1 = l/(fc). 

For a finite system, even though a true threshold does 
not exist, the crossover rewiring rate X C (N,L) increases 
as we increase the rewiring rate. For homogeneous net- 
works such as ER networks rewiring has little effect on 
the behavior of the disease since (k) 2 w (k 2 ). For het- 
erogeneous networks such as SF networks, the change is 
more significant. For SF networks with 7 > 3 we have 
argued that in the thermodynamic limit the threshold 
will increase continuously with the rewiring rate. On 
the other hand, for SF networks with 7 < 3 there is no 
threshold in the thermodynamic limit, except for an in- 
finite rewiring rate. 
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